Superfluidity and magnetism in multicomponent ultracold fermions 
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We study the interplay between superfluidity and magnetism in a multicomponent gas of ultracold 
fermions. Ward-Takahashi identities constrain possible mean-field states describing order parameters 
for both pairing and magnetization. The structure of global phase diagrams arises from competition 
among these states as functions of anisotropics in chemical potential, density, or interactions. They 
exhibit first and second order phase transition as well as multicritical points, met ast ability regions, 
and phase separation. We comment on experimental signatures in ultracold atoms. 
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Fermionic 5-wave superfluidity requires pairing be- 
tween fermions with different internal states. However, 
the nature of superfluidity ioi N = 2 component systems 
and > 3 is fundamentally different. Superfluidity is 
suppressed by magnetization ioi N = 2 because there is 
one way to pair and not every particle can find a partner. 
This led to proposals for exoticpaired states with bro- 
ken translational symmetry [2| or gapless excitations 
In contrast, superfluidity drives magnetization for 
A/" > 3. Here condensation energy favors enhancing the 
population for paired components to different degrees de- 
pending on the density of states and interaction energy. 
This occurs by cannibalizing the populations for unpaired 
components. 

Ultracold atoms offer direct access to multicomponent 
fermionic superfluids. Observation of two component, 
equal population superfluids used Feshbach resonances 
to tune interactions [l, 0, 0, Hi • The population for each 
component is both tunable and essentially conserved, en- 
abling later experiments on superfluidity with imbalance 
[§, [lO]- Recently, scattering lengths and locations of 
overlapping Feshbach resonances between all three of the 
nearly degenerate lowest sublevels of ^Li were measured 
[Tlj . This suggests experiments with N > are within 
reach. 

In this Letter, we consider the interplay of superfluidity 
and magnetism within a mean-field theory of > 3 mul- 
ticomponent fermions, each individually conserved. We 
derive a Ward-Takahashi (WT) identity (Eq. [7j) which 
provides fundamental constraints on possible mean-field 
states describing both pairing A and magnetization M. 
We demonstrate this WT identity naturally leads to a 
specific form for the microscopic pairing wavefunction 
which we call diagonal pairing states (DPS) illustrated 
in Fig. [TJ In these states, gapless excitations always 
exist for N odd and are still possible for N even. 

Having classified the mean-field states, we derive global 
phase diagrams describing the competition among the 
DPS as shown in Figs. [2] and [3l We focus on phase 
diagrams tuned by anisotropics in chemical potential or 
density but the structure is the same for anisotropics in 



FIG. 1: (color online) Pairing (ellipses) two components (col- 
ors) in all possible ways generates the diagonal pairing states 
as shown for N — ?>. Only these states generically satisfy the 
constraints on the microscopic pairing wavefunction imposed 
by Ward-Takahashi identities derived in the text, linear com- 
binations do not. Paired (unpaired) components have gapped 
(gapless) quasiparticle excitations. 



interactions. For A/" > 3, AA"*" acts as an external field 
for M through MAA^. This coupling vanishes identi- 
cally on group theoretical grounds for A^ = 2 where M 
couples to A only at higher order. The structure of the 
DPS shows pairing always drives magnetization for A^ 
odd and generically does so for A^ even through the cou- 
pling MAA"*". First order transitions and corresponding 
metastability as well as phase separation regions sepa- 
rate different DPS while second order transitions separate 
DPS and the normal state. These transitions terminate 
at bicritical and multicritical points. 

Previous theoretical works have also addressed super- 
fluditiy with A/" > 3 components of fermions. This in- 
cludes analysis of mean-field states for A/" = 3 [HI, [13] 
and A^ = 4 [3] as well as phase diagrams for A/" = 3 
[15I . [lit . We focus on classification of allowed pairing 
states through general symmetry arguments and Ward- 
Takahashi identities. This complementary approach al- 
lows us to obtain generic and robust features of both the 
resulting states and phase diagrams as well as provid- 
ing a unified perspective on pairing in multicomponent 
fermionic systems. 

We consider the action S{H) = Sq ^ SintiJ^) where 
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describing a dilute gas of fermions with mass m and 
attractive contact interactions where ipa, ipa with a = 
1 . . . N are Grassman variables. The partition function 



is Z{H,T) = JV[tl;,tlj]eyip\-J^dTjd^^S{H,T) 
is the inverse temperature, D the dimensionality, and 
ks = h = 1. 

We consider general H and T to derive WT iden- 
tities, but the physical system with individually con- 
served components has Ha/3 = ^a/3-f6 = 

J (^^7 + Sa6^(3-f — (7 ^ ^) contributing /ia^a and 
^a(3^a^(3 to the actiou where = i^a'^a- Here jiia is 
a chemical potential while Xap describes interactions be- 
tween fermionic densities. 

The action is U{N) symmetric for 11^ = 11 and Aq,/3 = 
A. This group acts as 

(3 (3 

with Ua(3 a unitary matrix. Anisotropics in ji^ and X^^ 
explicitly break U{N) U{1)^ describing N conserved 
densities in the normal state. The superfluid state spon- 
taneously breaks U{N) U{1)^ /7(1)^-^ where P 
is the number of non-zero pairing amplitudes. 

When Eq. [2] is viewed as a field redefinition 
transformation, the coupling constants also transform 

Hap H'^p = ^^s'^(^l^l5^5(3'> ra/375 ^'a(3-i5 = 

Y.^vpa^oc^iUpv'^^jypaUl^Ul^ in ordcr to describe the 
same physical system. By definition, Z(yH^ T) = 
Z{H',T') is invariant. 

This invariance arises from a field redefinition and is 
not a physical symmetry. However, it still gives a WT 
identity (see Ref. [l8|) expressing Z{H,T) = Z{H' X) 
under an infinitesimal transformation to first order 
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for arbitrary a, Equivalent ly, the above equation con- 
strains the expectation value of a redundant operator 
(Tot associated with the field redefinition of Eq. [2] to 
be zero. This constraint is most useful in the presence 
of anisotropics and is trivially satisfied in the symmetric 
limit. Analyzing infinitesimal transformations to second 
order may yield information about the symmetric limit. 
We classify possible mean-field states by analyzing the 
above constraint within the saddle-point approximation 
discussed next. 

Introducing bosonic pairing /S.a(3 and magnetization 
Ma(3 fields decouples the fermionic interaction through 
a Hubbard- Stratonovich transformation Mean-field 
states are given by zero-momentum, r-independent ex- 
pectation values for the order parameters Map-, ^ap sat- 
isfying the saddle-point equations 



Map = TgM E Gap, ^ap = ^^'^^ ^ 
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where cOn = (2n + l)7r//3 are Matsubara frequencies, k 
are momenta, T is the temperature and Xap = \gA,ap — 
\gM (1 + decoupling interaction anisotropics in 

the pairing channel. Normal Gap {i^n, k) and anomalous 
Fap (ia;n,k) Green's functions satisfy 



(5) 



where (,ap = (^^n — /'^^)^ap in a matrix notation with 
suppressed indices. As matrices, M^ = M and H^ = H 
are Hermitian while A-^ = — A is skew-symmetric. 
Diagonalizing these matrices gives 
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where [x\ is the fioor function, M^ (A^) real (complex) 
eigenvalues, and ujuj = vjvj = Sij ort honor mal eigen- 
vectors. From Map ~ (ipai^p), the physical interpre- 
tation of is the linear combination of fermions giv- 
ing magnetization M^. Similarly, Aap ^ {tpa'4^p) shows 
V2i-i, give the two linear combinations of fermions 
paired with amplitude A^. 

From Eq. El the WT identity relating M and A 
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(7) 

only constrains the off-diagonal elements as a matrix 
equation. In the normal state Aap = implying Map = 
for a 7^ p. Physically, U{1)^ symmetry of N conserved 
densities prohibits mixing between different components. 

Superfiuidity spontaneously breaks factors of U{1) 
through off-diagonal elements of Aap. We now show the 
WT identity requires the eigenvectors to be of the form 



(8) 



where Si^a is di N x N permutation matrix with exactly 
one non-zero matrix element in each row and column. 
We denote these states the DPS. Consider the generic 
case when /j.a 7^ /i/3, g^^ap ^ 9^^^' "^^^ ^ = 3, the DPS 
give all possible mean-field states. States not of DPS 
form do not satisfy Eq. [71 For A/" > 3, the DPS give 
generic mean-field states. States not of DPS form can 
in principle satisfy Eq. [71 but off-diagonal elements are 
generally overdetermined by the saddle-point equations 
implying non-zero values are inconsistent. 

We now discuss some properties of DPS. They simplify 
Eq. [H by decoupling them to \_N/ 2j saddle-point equa- 
tions each involving only A^, M2i_i, and M2i. M is also 
diagonal and respects U{1)^ symmetry as in the normal 
state while A breaks P factors of U{1) where P is the 
number of non-zero A^. The microscopic pairing wave- 
function for the DPS is given by using a Bardeen-Cooper- 
Schrieffer (BCS) 5-wave state to pair two and only two 
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FIG. 2: (color online) N = 3 phase diagrams for T > T^^^ 
(left) and T < T^^^ (right) against anisotropics in chemical 
potential fia (top) and density (bottom). NS denotes 
the normal state and {af3) a state with a, (3 paired. Solid 
(dashed) lines are first (second) order transitions while dotted 
lines are boundaries for met ast ability regions. MSi denotes 
regions with i additional metastable paired states and PSi a 
phase-separated mixture of i paired states. 



components at a time in all possible ways. The WT iden- 
tity prohibits states given by linear combinations of the 
DPS. This gives a discrete set of mean- field states with 
N\/P\{N-2P)\2^ DPS for fixed TV and P giving a total 
of {i/V2)^HN{-i/V2) distinct states for fixed N with 
Hn{x) the Hermite polynomial. 

Quasiparticle excitations are gapped for paired com- 
ponents and are gapless for unpaired components which 
remain a Fermi liquid. There are at most [N/2\ non-zero 
implying one component is always gapless for N odd. 
For N even, a fully gapped quasiparticle spectrum occurs 
only when all are non-zero. Refs. [l2|, [13] obtained 
similar results for = 3. 

For simplicity, we construct phase diagrams in the 
weak-coupling or BCS regime near the superfluid transi- 
tion for small anisotropics but stress the preceding anal- 
ysis of the DPS is general. We minimize the free energy 



F =Tr 
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given by F = — TlogZ where we take jj^a = Ep + (^/ia, 

9^,al3 = + = gMniEp) + SMa/S with 

Ep the Fermi energy, n^Ep) the free fermion density 
and from here on, quantities with a hat are rescaled 



with respect to T^^^ ^ the critical temperature without 
anisotropics. Here a, 6, c, d are Ginzburg-Landau param- 



eters with a^b^d ^ F 



^(D-2)/2 



proportional to the density 



of states at Ef describing particle-hole symmetric con- 
tributions. In contrast, c ~ Eli^ ^^^^ for D = 2 and 



c ~ E^^ ^^^^ logFi? for D 7^ 2 is essentially proportional 
to the derivative of the density of states at and de- 
scribes particle-hole symmetry breaking contributions. 

U (N) symmetric terms are under Tr, the matrix trace. 
Notice the 6M AA^ term where AA''' acts as an external 
field to SM. On group theoretical grounds, this term 
is non- vanishing only for N > 3. Terms outside the 
trace explicitly break U{N) /7(1)^, including a term 
quadratic in A. The structure of global phase diagrams 
follows from these two terms. 

Pairing AA''' drives magnetization SM unless AA''' ex 1 
is particle-hole symmetric with 1 the identity matrix. In 
this case, a shift in chemical potentials described by M 
yields no gain in condensation energy. Only when all 
components pair and \Ai\ = \Aj \ does this occur. For N 
odd, one component is always unpaired so magnetization 
always develops. For N even, only when the N/2 inde- 
pendent equations determining A^ give \Ai\ = \Aj\ is it 
possible to have pairing without magnetization. 

A similar situation occurs for p-wave pairing in ^He 
or the organic superconductors. Unitary states describe 
pairing decoupled from magnetization while non-unitary 
states describe pairing coupled to magnetization [20[. 
Only a single constraint |d x d*| = on the d- vector de- 
scribing p-wave pairing is necessary for a unitary state. 
It is more difficult for the N/ 2 independent equations de- 
termining Ai to give I Ail = \Aj \ for the analog of unitary 
states in multicomponent 5-wave pairing. 

We now discuss the = 3 phase diagram in Fig. [2] 
with given interactions satisfying the generic condition 
^9^a(3 7^ ^^A 7(5- consider both fixed chemical po- 
tential jla and fixed particle density 

TLqi as 



Ha = n(£;F)xo,a + niXi^c, + ^2X2,, 



(10) 



where xq = [1,1,1], xi = [1,-1,0]/V2, X2 
[l,l,-2]/V6 and /ii. Hi parameterize anisotropics. 

First consider fixed ft^ (top Fig. [2]). When T > 
small anisotropy favors the normal state. Increasing 
anisotropy favors pairing by increasing the density of 
states for some components at the expense of others. 
This drives a second order transition into one of the 
three DPS. Tuning the direction of the anisotropy drives 
first order transitions between DPS when two of these 
states are degenerate along lines of enhanced symmetry. 
First order lines terminate at bicritical points from which 
metastability regions where an additional DPS is locally 
stable branch out. This is the expected behavior for 
quadratic symmetry breaking |2l|, • Bicritical points 
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FIG. 3: (color online) Representative = 4 phase diagrams 
for T > T^^^ (left) and T < T^^ (middle) against in 
chemical potential fi^ with legend (right) for the global min- 
imum (top) and first metastable state (bottom). (ai/3i,...) 
denotes a^, Pi paired and white regions NS the normal state 
(top) or absence of metastable states (bottom). First (sec- 
ond) order transitions separate different paired states (paired 
states from the normal state). 

terminate at the U (3) symmetric multicritical point when 
T = Below first order lines separate DPS 

and a region near small anisotropy where all three DPS 
are (meta)stable appears. 

The Maxwell construction gives fixed phase dia- 
grams (bottom Fig. [2|) from those for fixed fta- The com- 
pressibility tensor tva/S oc —d'^F/dfiadfi/s at fixed fta has 
positive eigenvalues except at boundaries between DPS. 
Here, AA''' and thus SM jump discontinuously. First 
order lines at fixed fia when two DPS are degenerate ex- 
pand into phase-separated mixtures of those two states 
at fixed Ua- For T < T^^^, the zero anisotropy point 
where all three DPS are degenerate expands into a phase- 
separated mixture of those three states at fixed n^. 

Now consider the = 4 phase diagram in Fig. [3] with 
fixed chemical potential given by 

=^Fyo,a + Ao cos((9)yi,c,+ ^^^^ 
Ao sin((9) cos(0)y2,a + Ao sin((9) sin(0)y3,c, 

where yo = [1,1,1,1], yi = [1, -1, 0, 0]/V2, y2 = 
[l,l,-2,0]/V6, ys = [l,l,l,-3]/yT2 and i9, the 
anisotropics. 

For > 3, condensation energy favors development 
of more pairing amplitudes as T is lowered. The normal 
state competes with the six DPS with only one pairing 
amplitude for T > (top right Fig. [3]). However, 

the three DPS with two pairing amplitudes and all com- 
ponents paired eventually dominate the phase diagram 
for T < (top middle Fig. Phase diagrams stih 

exhibit second order transitions from the normal state to 
DPS and first order transitions between DPS. 

We now comment on trapping effects and detection 
methods for applications to ultracold atoms. In the large 



particle number BCS regime, the local density approxi- 
mation accurately maps phase diagrams at fixed chemical 
potential to phase diagrams with trapping. DPS are dis- 
tinguished by both densities and pairing amplitudes for 
different components. St ate- selective imaging of density 
distributions [9|, [lOl and radio-frequency spectroscopy of 
the pairing gap [23| can be used to detect signatures of 
the various paired states. 

In summary, we have studied the essential role 
magnetism plays in superfluidity of multicomponent 
fermions. By analyzing constraints imposed by Ward- 
Takahashi identities, we classified the allowed mean-field 
pairing states with both magnetization and pairing order 
parameters and used them to construct global phase dia- 
grams. These phase diagrams have a rich structure with 
first and second order transitions meeting at multicriti- 
cal points as well as metastability and phase separated 
regions. We discussed applications to ultracold fermions. 
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